Fractional and integer quantum Hall effects in the zeroth Landau level in graphene 
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Experiments on the fractional quantized Hall effect in the zeroth Landau level of graphene have 
revealed some striking differences between filling factors in the ranges < \u\ < 1 and 1 < \u\ < 2. 
We argue that these differences can be largely understood as a consequence of the effects of terms 
in the Hamiltonian which break 5(7(2) valley symmetry, which we find to be important for \v\ < 1 
but negligible for \v\ > 1. The effective absence of valley anisotropy for \u\ > 1 means that states 
with odd numerator, such as \u\ —5/3 or 7/5 can accommodate charged excitations in the form of 
large radius valley skyrmions, which should have a low energy cost, and may be easily induced by 
coupling to impurities. The absence of observed quantum Hall states at these fractions is likely due 
to the effects of valley skyrmions. For \u\ < 1, the anisotropy terms favor phases in which electrons 
occupy states with opposite spins, similar to the antiferromagnetic state previously hypothesized to 
be the ground state at v — 0. The anisotropy and Zeeman energies suppress large-area skyrmions, 
so that quantized Hall states can be observable at a set of fractions similar to those in GaAs two- 
dimensional electron systems. In a perpendicular magnetic field B, the competition between the 
Coulomb energy, which varies as B 1 ^ 2 , and the Zeeman energy, which varies as B, can explain the 
observation of apparent phase transitions as a function of B for fixed v, as transitions between 
states with different degrees of spin polarization. In addition to an analysis of various fractional 
states from this point of view, and an examination of the effects of disorder on valley skyrmions, 
we present new experimental data for the energy gaps at integer fillings v — and v = — 1, as a 
function of magnetic field, and we examine the possibility that valley- skyrmions may account for 
the smaller energy gaps observed at v = — 1. 

PACS numbers: 



I. INTRODUCTION 

The Dirac spectrum and high four-fold valley and spin 
degeneracy of Landau levels (LLs) in graphene lead to in- 
teresting phenomena in the quantum Hall effect (QHE) 
regime 1 . In particular, Coulomb interactions lift LL 
degeneracies, giving rise to new integer quantum Hall 
states with broken spin and/or valley symmetrjPHIH. 
Some of these states were conjectu red to have counter- 
propagating spin-filtered edge states^ESE^l, as well as un- 
usual spectra of collective excitations^. 

Recently, the progress in sample fabrication led to 
the observation of fractional quantum Hall (FQH) 
states in high-mobility graphene samples (suspended and 
graphene on BN}^"**', m DO th transport measurements, 
and in measurements of the local incompressibility us- 
ing a scanning single-electron transistor (SET) technique. 
Recent experimental studjEHl revealed an unexpected se- 
quence of FQH states in the zeroth LL. Certain FQH 
states in the interval —2 < v < — 1 of filling factors 
were missing or very weak, including states at filling 
factors v = — | , — | , — This is quite surprising, be- 
cause fractional states with denominator 3 are generally 
expected to be strong. In the interval of filling factors 
— 1 < v < 0, the FQH sequence was essentially iden- 
tical to that in found in GaAs-based 2D electron sys- 
tem^. More recently, incompressibility measurements in 
the range — 1 < v < have revealed transitions, as a 
function of the applied perpendicular magnetic field, at 



a fixed filling factor, which have been interpreted as tran- 
sitions between states with different spin polarizations^. 
It is worth noting that the absence of a v — —5/3 plateau 
was first observed in experimeniP^. Also, new multi- 
component FQH states arising due to valley and spin 
degeneracy of LLs in graphene were considered theoreti- 
callySH. 

Here we propose some possible explanations for the 
puzzling observations made in Refs. (19120122], We ar- 
gue that the key differences between the first partially 
filled Landau sub-level (—2 < v < —1) and the second 
one (— 1 < v < 0) arise from the greater importance in 
the second regime of the valley anisotropy terms in the 
Hamiltonian, which break the SU(2) symmetry with re- 
spect to the valley index. 

In what follows, we will use both the carrier filling 
factor v = n<&o/B±, where n is the net carrier density, 
B± is the magnetic field perpendicular to the sheet and 
$o is the flux quantum, and the "Landau level filling 
factor" / = v + 2, to describe filling of the zeroth LL. 
(When there are no carriers present, v = 0, the LL is 
half full, and / = 2). We shall concentrate on filling 
factors — 2 < v < 0, as the filling factors < v < 2 are 
related to these by particle-hole symmetry. 

The rest of the paper is organized as follows. In Sec- 
tion [n] we introduce the Hamiltonian and discuss how 
particle-hole symmetry can be used to relate the energies 
of different states. In Section IIIII we review the known 
results about integer states at filling factors v = — 1 and 
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0, and argue that the excitations of these states are dif- 
ferent, which makes the energy gap at v = —1 smaller 
than the gap at v = 0. In Section |IV| we study various 
FQH states in the filling factor interval —2 < v < — 1. 
We argue that the states at odd-numerator fractions gen- 
erally have skyrmion excitations, while the excitations of 
the even-numerator states are quasi-particles and quasi- 
holes. In Section[V]we turn to FQH states at — 1 < v < 0. 
We discuss the valley and spin structure of the states, 
finding that at all filling factors in this interval except 
for v — —1/3 there are two or more competing states 
with different degree of valley and spin imbalance. The 
competition between Coulomb energy and Zeeman en- 
ergy results in multiple phase transitions that occur as a 
function of magnetic field. In Section [VI] we consider the 
effects of disorder on various fractional and integer states, 
and argue that disorder can weaken and even suppress 
the odd- numerator states at filling factors — 2 < v < — 1. 
Section |VII| contains new experimental data on the en- 
ergy gaps of the v = — 1 and v — states, as well as 
a discussion about the nature of the excitations of these 
states. We summarize our findings in Secti on |VIII| and 
compare them to recent experiment d 19 * 20 * 2 ^. Finally, in 



the Zeeman interaction breaks spin rotation symmetry: 



Section IX we mention possible future directions for stud- 



ies of fractional quantum Hall effect in graphene and re- 
lated materials. 



II. PROJECTED HAMILTONIAN 

To understand the underlying symmetries, let us start 
with the microscopic Hamiltonian of the zeroth LL. We 
will neglect LL mixing, and consider projections of all 
operators onto the zeroth LL. For the analysis of spin and 
valley effects, it is convenient to separate the Hamiltonian 
into the leading 5J7(4)-symmetric part H sym and terms 
that break spin and valley symmetry. The symmetric 
part may be written as 
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sym 2 



- / drdv' : ft(r)ftr)V(r-r')ft(r')ftr') 



(1) 

where the symmetric interaction is given, at least approx- 
imately, by a Coulomb potential, V(r — r') = e | r e _ r ,| , 

where e is an effective dielectric constant. Also, ft (r) 
here is an operator that creates an electron in the zeroth 
LL. Since there are four fermion species, we view it as a 
four- component field, 

ft(r) = (^r),^O0.^W r )'^W r )) ( 2 ) 



H Z = E Z dvft{v)a z ftv) 



E z = 9 -p, B B 



(3) 



The most important valley-anisot ropic t erms arise due to 
lattice-scale Coulomb intera ction d 7 * 13 * 16 -*, as well as due to 
electron-phonon interaction d 14 * 15 [ There are two types of 
valley-anisotropic terms, which can be written as follows, 



dr 



ft(v)T z ftr) 
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(4) 



(5) 



The coupling constant g z originates from the electron- 
electron interactions, while g± includes contributions 
fro m bot h electron-electron and electron-phonon interac- 
tion d 7 * 14 * ^. The valley-anisotropic t erms c an get renor- 
malized and enhanced by LL mixing^Ell but here we 
will treat them as phenomenological parameters. We also 
assume that the interactions can be treated as instanta- 
neous in time. We note that LL mixing will also lead 
to some renormalization of the SU(4) invariant Coulomb 
interaction, on the scale of the magnetic length, but we 
ignore such effects here. 

The valley anisotropy terms in our model, Q and jHJ, 
have effect only when two electrons coincide in space. 
This reflects the fact that on the microscopic scale, sen- 
sivity to the electron valley or sublattice occurs only 
when two electron are separated by a distance of the 
order of the graphene lattice constant, which is much 
smaller than the magnetic length. We also remark that 
in the zeroth Landau level, the electrons in opposite val- 
leys are confined to opposite sublattices. 

With these assumptions the contributions of g z and 
g± to the energy per flux quantum, for a fixed value of 
v, will scale as g z /(2irl B ) and g±/(2Ttl B ), respectively. 
Here l B = y/hc/eB± is the magnetic length. (We will 
mostly consider the case when magnetic field is perpen- 
dicular to the sample, therefore, in the rest of the paper 
we identify B with B±). Consequently, this energy will 
scale the same way as the Zeeman energy, proportional 
to B, as previously remarked. (We note that the rela- 
tive importance of these terms can be altered, however, 
by application of a parallel field, which will increase the 
Zeeman energy but not affect the valley anisotropy ener- 
gies.) Because the energies due to Coulomb interactions 
scale as B 1 ' 2 , we see that the ratio between these and the 
valley anisotropy and Zeeman terms will change when the 
magnetic field is varied at fixed filling factor, and transi- 
tions between different ground state configurations could 
occur as a result. 



The four components may be regarded as forming a 
spinor in a four-dimensional space, which we shall refer to 
as "hyperspin space," and we shall refer to an arbitrary 
vector in this space as a "hyperspin" state. 

There are three main symmetry-breaking terms. First, 



A. Implications of electron-hole symmetry 

As was mentioned in the Introduction, the Hamilto- 
nian described above, when projected onto a single Lan- 
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dau level, has an exact particle hole symmetry. Specif- 
ically, it can be shown that given any many-body state 
|$) with LL filling factor / and energy E§, there exists 
a particle- hole conjugate state \^f), with LL filling factor 
4 — /, whose energy is given by 

/ — 2 

= E$ H — £?4 , (6) 

where E4 is the energy of the completely full Landau 
level, i.e., the state with / = 4. The above expression 
can be obtained by straightforward generalization of the 
analysis for the case of two Landau sub-levels performed 
in Ref. [23]. 

For the subsequent analysis, it is also helpful to note 
that there are additional electron-hole symmetries if we 
restrict the hyperspin states occupied by the electrons. 
For example, if we only consider many body states with 
electrons belonging to a single, fixed, hyperspin state, 
there is a particle-hole symmetry relating a state | $) with 
LL filling factor / and a state \^f) with LL filling factor 
1 — /, such that Ey = Eq, + (2/ — l)Ei. If we consider 
many-body states with electrons belonging to two speci- 
fied orthogonal hyperspin states, then there is a particle- 
hole symmetry relating a state |<&) with LL filling factor 
/ and a state |\&), with LL filling factor 2 — /, such that 

E 9 =E* + {f-\)E i . (7) 

We will use these relations below when discussing the 
energies of various fractional states. 

III. IQH STATES 
A. Nature of the ground state. 

Wc first consider the ground states at / = 1 and / = 2, 
which is to say at v = — 1 and v = 0. If the terms in the 
Hamiltonian which break the SU(4) symmetry are set 
equal to zero, the ground-state at / = 1 is predicted 
to be simple Slater determinant state, in which we have 
selected an arbitrary state \x) in the four-dimensional 
hyperspin vector space, and filled all orbital states in 
the Landau level with the given hyperspin. Similarly, at 
/ = 2, we expect the ground state to consist of two filled 
Landau levels with arbitrary orthogonal hyperspin states 
|Xi) and IX2) ■ (More precisely, it can be shown that these 
Slater determinants are exact ground states of the SU(4)- 
symmetric Hamiltonian for a wide range of interaction 
potentials, which almost certainly includes the standard 
Coulomb interaction. We note, however, that there could 
be lower energy states for more peculiar interactions.) 
The energies of these states are given precisely by the 
Hartree-Fock approximation. 

As long as the Hamiltonian is chosen to be SU(4) in- 
variant, the ground state energies of the Slater deter- 
minant states will be independent of the choice of the 
occupied hyperspin states. The valley anisotropy terms 



and the Zeeman energy lift the ground state degeneracy, 
at least partially. To find the true ground state or ground 
states, one should calculate the expectation value of the 
anisotropy terms for arbitrary choice of the occupied hy- 
perspin states, and choose a state which minimizes the 
total energy. 

For the case of / = 1, the energy minimization is quite 
simple. When only a single hyperspin state is occupied, 
there will be zero probability density for two electrons to 
sit at the same point in space because of the Pauli exclu- 
sion principle. Since the valley anisotropy terms in our 
Hamiltonian act non-trivially only if two electrons occupy 
the same point in space, their expectation value will be 
zero in any wave function that contains only a single oc- 
cupied hyperspin state . Therefore, splitting of the SU(4) 
degeneracy can only arise from the Zeeman term. The 
ground state will then have electron spin aligned with the 
applied magnetic field, but the valley pseudospin state \v) 
may point anywhere on the Bloch sphere. 

For / = 2, the valley anisotropy terms come into play, 
because the Pauli principle does not prevent electrons in 
different hyperspin states from coinciding in space. In 
this case ground states with several different symmetries 
are possible, depending on the signs and relative sizes of 
g z , g±, and the Zeeman energy. These states have been 
classified and their mean-field phase diagram is summa- 
rized in Ref. |16j . If the Zeeman term were large com- 
pared to the valley anisotropy terms, electrons in the 
lowest LL at v — would have their spins completely 
aligned with the magnetic field, and the ground state 
would thus be a valley singlet. However, if this were 
the case, the state at v = should have a finite electri- 
cal conductance, due t o the contribution of zero-energy 
energy ed ge stat e d 6 * 10 *, which is contrary to experimen- 
tal resultiPE^U. One possible state, which would have 
an energy g ap at the edge as well as in the bulk is the 
Kekule statePEH, which is a spin-singlet state in which 
all electrons occupy a single linear combination of the two 
valley states, with 50 % probability to be in either valley. 
Another possibility, which, however, seems unlikely in 
the presence of strong short-range repulsion, would be a 
"charge-denstiy-wave" state where electrons occupy both 
spin states in a single valley. This state was discussed in 
Refs. [71T3] . 

Another possibility, explored by Jung and MacDon- 
aloP^, as well as by Kharitonov-^, is an antiferromagnetic 
(AF) state in which electrons have spin in a direction si 
on one sublattice (say sublattice A) and have a different 
orientation s" 2 on the second sublattice. (Recall that for 
electrons in the lowest LL, electrons in valleys K and K' 
will reside on different graphene sublattices.) In the ab- 
sence of a Zeeman field, the orientations si and £2 will 
point opposite to each other, but can lie along an arbi- 
trary axis in space. In the presence of a weak Zeeman 
field, the spins will cant slightly and line up so that their 
sum si + §2 points in the direction of the magnetic field 
B, while Si — §2 lies in the plane perpendicular to B. 
As noted by Jung and MacdonalcP^ and subsequently 
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by Kharitono\ff^, this canted antiferromagnetic (CAF) 
state could transform continuously to a spin-aligned state 
in the presence of a sufficiently large parallel magnetic 
field. The AF state will be the ground state, in the ab- 
sence of a Zeeman term, if g z > 0, g± < 0, and g z > \g_i_\- 
The expectation value of the anisotropy energy in the AF 
state, per flux quantum, is equal to £"| n = —g z /(2wl 2 3 ). 
For small Zeeman coupling, the canting angle will be 
small, and the energy gain due to canting will be equal 

tO TTl 2 B E 2 z /\ gi ]. _ 

Current experimental evidence points to the CAF state 
being the ground state at v = in the absence of parallel 
field componenlPZl. At certain large parallel field, the in- 
sulating CAF state gives way to the metallic FM state^, 
similar to the case of bilayer graphene 2 -. To characterize 
the relative strength of the valley anisotropy and Zeeman 
interactions, we introduce a quantity 

(8) 

\g±\ 

Based on experimental evidence, including the behavior 
of transport coefficients at v = in a tilted magnetic 
field, it appears that 7 is of order 0.1^. Therefore, the 
Zeeman energy is small compared to the valley anisotropy 
energy at / = 2, and the canting angle is small. 

In our analysis of integer and fractional states below, 
we will consider parameter values that favor CAF state 
(<? z > 0,<?j_ < 0, and g z > \g±\), and also mostly assume 
that the Zeeman interaction is much smaller than the 
valley anisotropy. 

B. Charged excitations at / = 1. 

Given our characterizations of the ground states at / = 
1 and / = 2, it is natural to ask about the elementary 
charged excitations in the two cases. As we shall see 
below, the nature of excitations is quite different in these 
two cases: at / = 1, the excitations are skyrmions, while 
at / = 2 the excitations are electrons and holes. 

We first consider the case of / = 1. The simplest model 
for a hole excitation is the state in which one simply 
removes one electron from the filled level with the chosen 
hyperspin state |x). The simplest model for an electron 
excitation would be to add a single electron in one of the 
hyperspin states orthogonal to |x). The energy of a hole 
excitation is given by 

A . pir e 2 

irrespective of its hyperspin. Quantity Ao is the energy 
gap for creating an electron-hole pair in the absence of 
symmetry-breaking terms. In the case of graphene, the 
energy of an electron excitation, however, contains a con- 
tribution from valley-anisotropic terms in the Hamilto- 
nian, and therefore is sensitive to the hyperspin of the 
extra electron. Let us assume that the ground state at 



/ = 1 is polarized along K f direction in the hyperspin 
space. Then we can consider extra electrons with hyper- 
spin K' l or K' f. Their respective energies are given 
by 

Ee ' AF ~T~2^| + ^' Ee ' FM ~2 2^Z| Z ' 

We also note that adding extra electron with K I will 
cost more energy. Given that ^> Ez, the electron 

excitation with hyperspin K' I (" AF" one) has the lowest 
energy. 

For Coulomb interactions, however ,we expect that the 
lowest energy charged excitations at / = 1 will be 
skyrmions, in which the hyp erspi n is misaligned rela- 
tive to |x) over a large are a 25 ! 29 ! (for a discussion of 
skyrmions in SU(N)-symmetric quantum Hall ferromag- 
net, see Ref. [30]). In the absence of symmetry breaking 
terms, Coulomb interactions will favor skyrmions of in- 
finite radius, where the skyrmion energy is found to be 
one-half of the energy of a single particle electron or hole 
excitatior l 25 l 29 [ In GaAs, the radius of the area of the 
skyrmion, and the number of overturned spins, is lim- 
ited by the Zeeman energy, which gives an energy cost 
proportional to the skyrmion area^l. 

In graphene, hyperspin can take 4 different values, so 
several kinds of skyrmions are possible. One possibility 
is to consider valley skyrmions - textures in which valley 
polarization is smoothy varied but spin is aligned with 
the magnetic field. There is no Zeeman cost for such a 
valley skyrmion, and we would expect that the size of an 
isolated skyrmion in an ideal sample would be arbitrarily 
large, to minimize the Coulomb interactions. (In prac- 
tice, the size of a skyrmion might be limited by disorder 
or by the distance between skyrmions). In any case, the 
energy of a valley skyrmion should be close to 1 /2 of that 
of a single electron or hole excitation: 

E sk = ^±E Z + E V , (9) 

where E v denotes the contribution of valley-anisotropic 
terms, and ± refers to hole (electron) skyrmions. 

To examine the effect of valley anisotropy terms, first 
note that a large area skymion may be characterized by 
a space dependent orientation in hyperspin space for the 
occupied electron state. At the center of the skyrmion, 
the hyperspin should point in a direction determined by 
a spinor \x2) that is orthogonal to the state \x) which 
characterizes the background far from the skyrmion. We 
define a normalized probability density P for two elec- 
trons to coincide at a point r, by 

P(r) = (27r/|) 2 (: ft(T)fa)ft{T)j>(T) :) (10) 

Because the valley anisotropy terms are relevant only 
when two electrons coincide in space, any contribution of 
these terms to the energy density of the skyrmion must 
be proportional to P(r). (We note that two electrons 
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with valley indices K, K' in the zeroth LL can never re- 
side on the same lattice site, because the zeroth LL wave 
functions in the K(K') valleys reside solely on the A(B) 
sublattice. So, by "coincide in space" we mean that two 
electrons are much closer than the magnetic length lg.) 
A priori, we would expect that P will be of order l B / R 2 
inside the skyrmion of radius R (one extra electron spread 
out over area R 2 ) and that it will fall to zero outside. 

A deeper analysis of the skyrmion structure at / = 1 
leads to the conclusion that P(r) = in the case of the 
hole skyrmiorPD at least in the limit of large skyrmion 
radius. Therefore there is no anisotropy energy for the 
hole skyrmion, and 

Ao 
4 

For the electron skyrmion, P^O, since we have added an 
extra electron. Furthermore, it follows from particle-hole 
symmetry that J P(r) d 2 r = 1 for the electron skyrmion. 
Therefore, there will be a non-zero anisotropy energy in 
this case, equal to E v = — 29 „ ± t A 9z . This gives the energy 
of the spin-aligned electron valley skyrmion: 



P. 



k,h 



Ez- 



V FM _ 
^sk,e — 



Ao 
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-E z - 



2ffi 



fJz 



2nl 2 B 



(11) 



This is the energy in the limit of infinite skyrmion 
size. Smaller skyrmions will be more costly due to the 
Coulomb interactions. 

As an alternative to the valley skyrmion, one may con- 
sider a spin-involved skyrmion, where \2 bas a spin ori- 
entation opposite to the magnetic field. In this case, \x2) 
may have a valley pseudospin that is the same or is differ- 
ent from that of |x)- Taking into account the implications 



of Eq.fl7 
is the 



, we find that if the lowest energy state at / = 2 
state, the electron skyrmion at / = 1 with low- 
est anisotropy energy will be one where the hyperspin 
state far from the skyrmion has electrons in one valley 
and spin aligned with the magnetic field, while at the 
center of the skyrmion the electrons sit in the opposite 
valley and have reversed spin. The anisotropy energy in 
this case would be E^ m 2 = — (2irl%)^ 1 (g z ), just like 
for the single electron-like excitation. 

The total energy of the "AF" skyrmion is given by: 

Ao 
4 



pAF _ 
^sk,e — 



2nl 2 B 



(2n s - 1)E Z + D 



(12) 



where R is the nominal skyrmion radius, the term 
De 2 /eR arises from the long-range Coulomb interaction 
energy (D is a numerical constant of order unity), and 
n s (R), the number of electrons with reversed spins, is 
proportional to (R/Ib) 2 , multiplied by a constant that 
depends logarithmically on R and Ez- The Zeeman inter- 
action favors small skyrmions, while the Coulomb inter- 
action favors infinitely large skyrmions. The competition 
of these two terms le ads to an optimal intermediate-size 
skyrmion with energjPsM] 

Ao _ g z 
4 2tt1 2 b 



rpAF 
^sk,e 



^A ^ /3 |log(10gz)r/ 3 . (13) 



Here 



and 



3 4/3 5/6 



E? 



gz 



e 2 /eh 



is the Zeeman interaction in units of Coulomb interaction 
energy (notice that the dimensionless Zeeman interaction 
differs from the definition adopted in Refs. 29 32] by a 
factor of 4). The expression (131 is valid in the limit 
gz <C 1. The number of electrons with reversed spins 
may be found from 



(2n s - 1) 



dE z 



~g~ 2/3 \\og(l0~g z )\ 



1/3 



The difference of the energies of two kinds of skyrmions 
is given by: 



AP, 



pAF _ pFM 
sk,e sk.e 



AA .^ /3 |log(10£z)| 



1/3. 



%±l 

2ttI% 



(14) 

The first term scales as P 2 / 3 with magnetic field, while 
the second term scales as P. Therefore, at low fields we 
expect the large valley skyrmions to be the lowest-energy 
excitations, while at high fields, spin-reversed skyrmions 
will have lower energy. The critical field at which the na- 
ture of the skyrmions changes is found from the condition 
APg/t.e = 0, which can be rewritten as follows, 



A^\~9 Z \-^\iog~ 9z \y 3 



This gives the critical dimensionless ^-factor (with loga- 
rithmic precision): 



9z 



A'y fir 
~Y\2 



3/2 



log 1 / 2 



10 



3/2' 



For 7 = 0.1, and e = 5.24 (the RPA result for intrinsic 
screening function of graphene^) , this yields 

gz » 0.0075, P* « 1.8 T. 

Thus, at P > P* we expect the lowest-energy excitations 
to be spin-reversed skyrmions. 

The number of reversed spins for AF skyrmions can 
be accurately determined using numerical results of 
Ref. |34I35| (the effective sigma-model result is valid only 
in the limit when parameter gz is very small, which is 
not necessarily the case for experimentally relevant fields 
P ~ 10 T). Within a classical (i.e., Hartree-Fock) de- 
scription, the mean number of electrons with reversed 
spins n s in a skyrmion will vary continuously with the 
radius of the skyrmion. In a proper quantum descrip- 
tion, however, the skyrmion must be an eigenstate of the 
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spin angular momentum parallel to the applied magnetic 
field, and hence the value of n s will be quantized and have 
only integer values. Similarly, the number of electrons n v 
with the reversed valley index K must also be an integer. 
Different types of skyrmions are thus properly character- 
ized by specifying the two integers n s and n v . For the AF 
skyrmion described above, we have n v = n s , while for the 
pure valley skyrmion, we have n s = 0, n v = oo. We may 
also consider skyrmions in the quantum equivalent of a 
canted AF state, which would have < n s < n v < oo. 
However, it is not clear that such skyrmions would be 
favored in the parameter range of interest. 

For magnetic fields B = 2.5,5,10T we find, using re- 
sults of Ref . [35] . n s « 5,4,3 for AF skyrmions. The 
number of reversed spins can be measured in parallel field 
experiments. Recent transport experiments on samples 
on BN substrate found that the activation gaps at / = 1 
were sensitive to the parallel field component, indicating 
that they involve spin flipa^. The effective ^-factor was 
found to be somewhat enhanced, which is consistent with 
small spin skyrmion excitations. It should be noted, how- 
ever, that the spin-reversed and valley skyrmions might 
have different sensitivity to disorder, and their trans- 
port properties might be different as well. Therefore it 
is possible that large valley skyrmions are present and 
contribute to compressibility, but are more easily pinned 
by disorder and therefore their contribution to transport 
is suppressed. Future compressibility measurements in 
tilted field might shed light on the nature of excitations 
at / = 1. 



C. Charged excitations at / = 2. 



imize Zeeman interaction) has an energy 



E K — An — Ez — 



2g± 

2tW| 



The corresponding activation gap is then given by 

A„ =0 = E h +E e = 2A - 2E Z + 9z ~l 9± . (15) 

Turning to the possibility of having skyrmions at v — 
0, we note that several types of charged skyrmions are 
possible. A charged skyrmion might have a filled LL for 
K j" and a valley-pseudospin texture which mixes K' J, 
and K \.. At the center of the skyrmion the system would 
be locally in a CDW state, with both spin states in valley 
A, which would cost an energy proportional to g z and to 
the area of the skyrmion. Similarly, we could consider a 
texture that mixes K' \. and K' f. Such a configuration 
would cost an energy proportional to \g±\, multiplied by 
the skyrmion area. Therefore, the energy of the charged 
skyrmion involves a valley anisotropy cost proportional 
to the skyrmion area. 

Given that the valley anisotropy is large, we there- 
fore expect the skyrmions at / = 2 to be unfavorable. 
The relevant parameter that determines the size of the 
skyrmions and their energy is g — ^^tj^ 3 ^- Tak- 
ing 7 = 0.1, we obtain g w 0.35, which is well above 
the critical values g ~ 0.06 at which skyrmions become 
completely su ppress ed and lowest energy excitations are 
electrons/holesPESl. Thus we conclude that the lowest- 
energy excitations at / = are single holes and electrons, 
with the corresponding gap given by Eq.( 15 1. 



In an SU(4)-symmetric model, skyrmions are possible 
at filling fraction / = 2, and have the same energy as at 
/ = 1. In the case of graphene, however, the short-range 
valley-anisotropic terms presumably favor an AF state, 
and, in the presence of Zeeman energy, the symmetry is 
broken down to Z(2) x {7(1) (where Z{2) is the valley 
symmetry and U(l) is the symmetry related to the rota- 
tion of the spin orientation in the xy-plane). Because of 
the low symmetry, skyrmions cost a large anisotropy en- 
ergy proportional to their area, in contrast to the case of 
/ = 1. As a result, skyrmions are not favorable and the 
lowest excitations are single-particle- like. We expect that 
this would make the energy gaps at / = 2 significantly 
larger than the gaps at / = 1. 

Let us suppose that the ground state at v — is the 
AF state, with occupied states K f and K' L A single 
hole excitation, in which we remove an electron with spin 
that is anti-aligned with magnetic field (hyperspin K' J,) , 
has an energy 

E h = A -E z + ^, 
while an electron excitation with hyperspin K' f (to min- 



IV. FQH STATES FOR < / < 1 

We start our analysis of the FQH effect in graphene by 
considering various states in the interval of filling frac- 
tions < / < 1. 

A. / = l/3 (u=-S/3) 

The ground state at / = 1/3 is presumably well de- 
scribed by the Laughlin trial wave function, with the 
electrons in a single hyperspin state, with spin aligned 
by the magnetic field and arbitrary orientation of the val- 
ley pseudospin. The Laughlin state is the exact ground 
state in the limit of a zero-range interaction potential, 
and is known to be very accurate for a Coulomb poten- 
tial. Moreover, the assumption that the ground state 
contains only one hyperspin state is consistent, as the 
valley-anisotropy terms have no effect on such a state 
due to the Pauli exclusion principle, and the state has 
been already chosen to minimize the Zeeman energy. 

The lowest energy charged excitations at / = 1/3 are 
probably valley skyrmions. The size of such a skyrmion 
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can be very large, as it does not cost any Zeeman en- 
ergy, in contrast with the case of a conventional semi- 
conductor, where the skyrmion involves spin orientations 
that are not aligned with the magnetic field. As far as 
we arc aware, there arc currently no reliable calculations 
of the spin-stiffness constant for an / = 1/3 state with 
Coulomb interactions, so it is difficult to estimate the 
size of the energy gap that would be produced. H oweve r, 
based on numerical calculations for small systems^HZl it 
is expected that the skyrmion energy will be significantly 
smaller than the energy of a Laughlin quasiparticle or 
quasihole. 

Unlike the skyrmions at / = 1, the excitations at 
/ = 1/3 should have only a very small probability of 
two electrons being at the same position in space, for 
quasi-electron-like as well as a quasi-hole-like skyrmions. 
Indeed, it is possible to construct a trial wave function for 
a pseudospin-reversed electron, or for a skyrmion, which 
costs energy due to the V% pseudopotential, but has no 
contribution from Vq, and it is probable that even for 
Coulomb interactions, the skyrmion wave func tion will 
tend to strongly avoid such a contributio n 1 38 ! 39 !. Conse- 
quently, the valley-anisotropic interactions should have 
very little effect on the energy of the skyrmions. 



B. / = 2/5 (u = -8/5) 

We expect that the ground state at / = 2/5 should be 
a pseudospin singlet state, with two valley states equally 
occupied and with electron spins aligned with the mag- 
netic field. The elementary trial wave function for this 
state is the Halperin 332 wave function^, in which the 
amplitude vanishes as the cube of the separation when 
two electrons in the same valley state approach each 
other, and vanishes quadratically for electrons in oppo- 
site valleys. This is an exact ground state for a system 
of N > 2 species of particles with an SU(N)-invariant 
Hamiltonian in the limit of a zero-range repulsive inter- 
action, since it is the only wave function at filling 2/5 in 
the lowest LL which as no contributions from the Vq and 
Vi pseudopotentials. It is also found to be an excellent 
approximation in the case of a Coulomb potential. 

In GaAs, the spin-singlet state is in competition with a 
spin-aligned 2/5 state that has a higher Coulomb energy 
but is favored by the Zeeman term, so there exists a tran- 
sition between the two states as one varies the magnetic 
field at fixed filling factor, or if one increases the Zee- 
man energy by application of a parallel magnetic field. 
In graphene, however, the spins are already fully aligned 
with the magnetic field in the valley-singlet state, so there 
is no Zeeman penalty for forming the state. Therefore, 
we predict that a state where all electrons have the same 
hyperspin configuration should never be the ground state 
in graphene at / = 2/5, and there should be no transition 
as one varies the magnetic field. 



C. / = 2/3 (i/ = -4/3) 

At first glance, there are two plausible candidates for 
the state at / = 2/3. One may consider a state in which 
all electrons have the same hyperspin configuration, or 
one may consider a state which is a valley- singlet but 
spins aligned by the magnetic fielcP^, as for the / = 2/5 
state discussed above. In GaAs, at / = 2/3, it is found 
that the spin singlet state has a lower Coulomb energy 
than the spin-aligned state, but as at / = 2/5, the spin- 
aligned state is favored by the Zeeman energy, and there 
is a transition between the two states depending on the 
ratio of the Zeeman energy to the Coulomb energy energy 
scale. In graphene, by contrast, the valley singlet and the 
hyperspin-aligned states would have the same Zeeman 
energy, so we would expect the valley singlet to be the 
ground state at all magnetic fields. However, we should 
also consider the possible contribution of the anisotropy 
terms before concluding that this is the case. 

As noted previously, a state with all electrons in the 
same hyperspin state will have no energy contribution 
arising from the valley anisotropy terms, due to the Pauli 
exclusion principle. On the other hand, the valley singlet 
state at f — 2/3 could have a non-zero probability den- 
sity P for two electrons in opposite valleys to coincide at 
the same point in space. [We continue to use the normal- 
ization of P defined by ( fio| ), which is such that P = 1 in 
the AF state at f = 2] Then if P 2 /3 is t ne value of P in 
the / = 2/3 valley-singlet state, the anisotropy energy, 
per flux quantum will be — {g z + 25'x)-F > 2/3/(27rZ^), which 
will be positive if (g z + 2g±) < 0. This inequality is con- 
sistent with the condition g z > —g± > for stability of 
the AF state at v = 0, but is not required by it. As pre- 
viously noted, the anisotropy energy per electron at fixed 
/ will scale linearly with the applied magnetic field, just 
as for the Zeeman energy in GaAs, so if (g z +2g±) > 0, it 
is conceivable, in principle, that there could be a transi- 
tion between the two types of ground states in graphene 
at / = 2/3. However, it appears that the value of P 2 /3 
is actually extremely small, so that the valley singlet is 
always the ground state. 

Unlike the case of the 332 wave function at / = 2/5, 
there is no simple trial wave function for the valley-singlct 
state at / = 2/3 where one can evaluate the value of 
P. Nevertheless, numerical calculations on a finite sys- 
tem suggest that P 2 /3 nas a value less than 10 -4 for 
pure Coulomb interactional. This is, in fact, very much 
smaller than the value P = 1/9 that one would obtain 
if one filled one third of the available states in each val- 
ley with spin parallel to the magnetic field, and assumed 
that there were no spatial correlations between electrons 
in opposite valleys. The small value of P2/3 also guaran- 
tees that any advantage in anisotropy energy could never 
outweigh the Zeeman penalty for a state in which spins 
are not all aligned with the magnetic field. If we add 
a skyrmion to the valley singlet state, it will necessarily 
cost a substantial Zeeman energy, so we do not expect 
there to be low-energy charged excitations in the form of 



large-radius skyrmions. 

We note that if the ground state at / = 2/3 were actu- 
ally the state with electrons confined to a single hyperspin 
state, it would be a particle hole conjugate of the state 
at / = 1/3. Just as for the / = 1/3 state, there would be 
low-energy skyrmion excitations, which mix in the second 
valley state but do not cost any Zeeman energy. 

Within the language of composite fermions, the valley- 
singlet state at / = 2/3 would be described as a state 
where one has doubly filled the lowest composite fermion 
LL, in an effective magnetic field that is opposite in sign 
to the applied field, with particles in two orthogonal hy- 
perspin states. Alternatively, it is possible to describe 
the valley-singlet 2/3 state as one that can be obtained 
by adding a 1/3 density of hole-like charged excitations 
to the fully polarized state at / = 1. In this case, the 
charged excitations would be minimally-sized hole-like 
skyrmions, which may be described as a bound state of 
two missing electrons in the occupied hyperspin state and 
one added electron in the reversed valley state. As the 
difference in occupation numbers of the two valley states 
is decreased by three for each added hole, the difference 
will be reduced to zero at / = 2/3. 



Laughlin m = 5 state formed from hole-type skyrmions 
of the type described above. 

This state will be fully spin-aligned but should be in- 
variant under SU(2) rotations in valley space. There- 
fore, it cannot accommodate pure valley skyrmions. Any 
skyrmions must therefore involve the spin orientation, 
with an energy cost due to the Zeeman field. Conse- 
quently, the state could have a substantial energy gap 
for charged excitations. 

Similar arguments suggest that the ground state at 
/ = 6/7 should be a valley singlet state, and that it might 
also be a quantized Hall state. However, a quantized 
Hall state has not been observed experimentally at this 
fraction. 

Unpublished calculations of finite systems by R. MorP2l 
suggest that for pure Coulomb interactions, spin-singlet 
ground states should exist at / = 4/5 and / = 6/7, 
and that the probability density P for two electrons to 
coincide is very small in both of these states, of order 
1CP 3 . This implies that the effects of valley anisotropy 
should be negligible in these states, and it is consistent 
with the skyrmion construction starting from / = 1 state, 
as outlined above. 



D. / = 4/5 and 6/7 (v = -6/5, -8/7) 

We may consder several possible candidates for a state 
at / = 4/5. One possibility would be a completely hyper- 
spin polarized state, in which all electrons occupied the 
same valley state and had spin aligned with the magnetic 
field. Such a state would be a particle-hole conjugate of 
the state at / = 1/5, and would be described as a Laugh- 
lin 1/5 state of holes in the full hyperspin-polarized Lan- 
dau level. This state would presumably be the correct 
ground state in a model of short-range interactions, with 
a dominant Vq repulsive pseudopotential, but it is un- 
likely to be the true ground state in the case of Coulomb 
interactions. 

As was previously noted, at / = 1, the lowest en- 
ergy hole excitation, in the case of Coulomb interactions, 
should be a large-area skyrmion with many reversed val- 
ley spins, rather than a single missing electron from the 
full polarized LL. Thus it seems likely that the ground 
state at / = 4/5 should be constructed out of skyrmions 
rather than of holes in the polarized LL. In particular, let 
us consider a finite size skyrmion consisting of four miss- 
ing electrons bound to three electrons in the reversed 
valley state. If we add such skyrmions to the / = 1 
state, at a density of one skyrmion for every five flux 
quanta, we obtain a state with / = 4/5 and equal popu- 
lation of the two valley states. One possibility is that the 
skyrmions would form a Wigner crystal, in which case 
the Hall conductance would lie on the v = — 1 plateau, 
and there would not be an FQH state at / = 4/5. How- 
ever, recent experiments do show clear evidence of an 
FQH state at / = 4/5. Therefore, we believe that the 
FQH state at / = 4/5 can be understood, roughly, as a 



E. Other fractions with < / < 1 

At / = 4/9, we see two candidate ground states. One 
possible state would have equal numbers of electrons in 
each of two hyperspin states, occupying states in both 
valleys with spins parallel to the magnetic field. The 
ground state in this case could be described in the com- 
posite fermion picture as a state with two types of com- 
posite fermions occupying the lowest two effective LLs. 

An alternative state could have electrons in all four 
hyperspin states. This state could be described by a 
generalization of the Halperin 332 state, in which the 
many-body wave function contains factors that vanish as 
(zi — Zj) 3 for particles in the same hyperspin state, but 
vanish as [zi — Zj) 2 for particles with different hyperspin 
states (here z% — Xi + iyi is the complex coordinate of 
ith particle). Such a state would presumably minimize 
the Coulomb energy but would pay a price in Zeeman 
energy. In principle, there could be a phase transition 
between the two types of states, if one varies the mag- 
netic field at fixed LL filling. An FQH state at / = 4/9 
has been observed in experiments, but no phase tran- 
sition has been observed as yet. We also note that one 
could consider other states at / = 4/9, e.g., ones in which 
three different hyperspin states are occupied. However, 
given the double degeneracy of CF levels at / = 4/9, it 
seems unlikely that such states would be favorable. 

In a similar way, at / = 3/7, we might consider a 
ground state with three different hyperspins occupied, 
each with 1/7 electron per flux quantum, or a state with 
all electrons in one hyperspin state, or a state with 2/7 
electron per flux quantum in one hyperspin state, and 
1/7 in another. Similarly, at / = 3/5, we could consider 
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states with three different hyperspins occupied, with all 
electrons in one hyperspin state, or having two hyper- 
spin states with filling factors of 2/5 and 1/5. All of 
these states would have charged skrymion excitations 
that cost no Zeeman energy, and little or no anisotropy 
energy. Therefore, we might expect them to have only 
a small energy gap, and the FQH state would be easily 
suppressed by disorder or finite temperature. A FQH 
states at / = 3/7 has not been observed yet, and there is 
only a faint hint of an incompressible state at / = 3/3^. 



V. FQH STATES FOR 1< / < 2 

In the discussion below, we shall only consider can- 
didate ground states in which only two fixed hyperspin 
states, \xi) and \x2) may be occupied. (The remaining 
two hyperspin states are completely empty.) States of 
this type are consistent with the symmetry of our Hamil- 
tonian, and it seems like a reasonable assumption that 
the actual ground states will satisfy this restriction, but 
we do not have any rigorous argument that this must be 
the case at all filling factors. 

Following our assumptions, the ground state at LL fill- 
ing / in the range 1 < / < 2 should be related by particle- 
hole conjugation to the ground state one would obtain for 
the same choice of |xi) and |%2) at a LL filling /' = 2 — /, 
in the range < /' < 1. In particular, the expectation 
value of the anisotropy energies in the states at / and 
2 — /, for the given choice of hyperspin states, would be 
related, following Eq. by 



/ 



(/ - w 



(16) 



where -Ef n is the anisotropy energy of the state at / = 2 
for this choice of |xi) and |%2)- However, the conjugate 
state at /' may not be the true ground state at that filling 
factor, because the optimal choice of hyperspin states 
may be different at /' than at /. We found previously, 
that the optimum ground states for < / < 1 were 
obtained when the two constituent hyperspin states had 
opposite valley indices but both spins aligned with the 
magnetic field. However, if the Zeeman energy is small 
compared to the anisotropy terms at / = 2, so that the 
AF state is the ground state at that filling factor, we may 
guess that for 1 < / < 2 we should choose the hyperspin 
states to be the same ones that are represented in the AF 
state, e.g., K f and K' \., so the system is not fully spin- 
polarized. (We assume, for the moment, that the Zeeman 
field is small compared to the anisotropy energy in the 
AF state, so we may ignore canting of the spins. We shall 
discuss effects of canting in a separate subsection below.) 

In the following discussion we shall assume that \xx) 
and |x2) are given by K f and K' |, for 1 < / < 2. With 
this choice, we see that i?| n is the anisotropy energy of 
the AF state at / = 2, given by E% n = -g z /(2nl 2 B ) < 0. 



A. / = 5/3 (y = -1/3) 

The obvious candidate for this filling fraction is a maxi- 
mally polarized state |<i>), where we align the spins quan- 
tization axis with the magnetic field, and remove from 
the / = 2 AF state one-third of the electrons in the hy- 
perspin state with spin opposite to the magnetic field 
(K' I), while retaining the completely full LL of elec- 
trons in the hyperspin state K f favored by the Zee- 
man interaction. We may describe this state as a Laugh- 
lin state of hyperspin-aligned holes in the / = 2 state. 
The particle-hole conjugate of this state would be an 
/ = 1/3 state |$') where all electrons are in the hy- 
perspin state K' J.. That state will have no anisotropy 
energy, because of the Pauli exclusion principle, so the 
proposed state at / = 5/3 should have an anisotropy 
energy E^ 3 = (2/3)£f n , according to (l6|). 

We argued, previously, that the true ground state at 
/ = 1/3 should have all electrons in a single hyperspin 
state with spin parallel to the Zeeman field (here, spin f) 
and arbitrary choice for the valley pseudospin state. If 
the Zeeman energy were absent, the state |<3>') would have 
the same energy as the true ground state at f — 1/3, so 
it must also be a state that minimizes the combination 
of Coulomb and valley-anisotropy energies. It follows 
by particle- hole symmetry, that at / = 5/3 the wave 
function |$) must also minimize these energies, among 
states within the class of we are considering, in which 
electrons are restricted to the hyperspin states K j" and 
K' \,. As the state |<i>) also minimizes the Zeeman energy, 
it is the correct ground state at / = 5/3, at least within 
the class we are considering. 

Although the ground state energy at / = 5/3 is directly 
related to the energy of the ground state at / = 1/3, the 
excitation energies will be quite different. In particular, 
a skyrmion at / — 5/3 will inevitably engender a cost in 
Zeeman energy and/or anisotropy energy. Consequently, 
we expect that there will be a substantial energy gap at 
/ = 5/3 and a well established quantized Hall state, in 
contrast with the situation at / = 1/3. 



B. / = 4/3 (i/ = -2/3) (/ = 4/3) 

If we restrict our considerations to many-body states 
with only two occupied hyperspin states, the ground state 
at / = 4/3 should be related by particle-hole conjugation 
to the ground state at / = 2/3 for the same choice of hy- 
perspin states. However, the optimal choice of the hyper- 
spin states should be different in the two cases. We found 
previously that the optimum ground state at / = 2/3 was 
obtained when the two constituent hyperspin states had 
opposite valley indices but both spins aligned with the 
magnetic field. However, at f — 4/3, if the Zeeman en- 
ergy is small compared to the anisotropy terms, and the 
AF state is the ground state at / = 2, we may guess that 
it would be better to choose the hyperspin states to be 
the ones that are represented in the AF state, that is K f 
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and K' J, (we will justify this below, in the subsection on 
spin canting). Now, if we choose the spin quantization 
axis to be parallel to the Zeeman field, the state with fully 
polarized holes can have a lower Zeeman energy than the 
state with equal populations of the two hyperspin states. 
Thus the Zeeman field can now drive a phase transition 
between the two states. 

To make this more precise, we must again consider 
the anisotropy energy, which will be non-zero in both of 
the states we consider. From Eq.Q, we see that the 
anisotropy energy of a state at / = 4/3 is related to that 
of the conjugate at / — 2/3 by 

EX/z = Ey 3 + (l/3)£ 2 an . (17) 

Since there can be no anisotropy energy in the fully po- 
larized state at / = 2/3, we see that anisotropy energy 
for the maximally polarized state at / = 4/3 for the same 
value of the magnetic field, will be equal to one-third of 
the anisotropy energy at / = 2. The anisotropy energy 
for an / = 2/3 state with equal occupations of K t and 
K' I would be given by E$J 3 = g z P 2 / 3 = P2/3E™, where 
P2/3 is the probability density for finding two electrons 
at the same position in space in the / = 2/3 state. The 
value of -P2/3 does not depend on the choice of hyperspin 
states, as the orbital part of the wave function is deter- 
mined by the SU(4) invariant part of the Hamiltonian, 
and we have previously argued that the value of P2/3 
should be very small. Consequently, we expect that the 
difference in the anisotropy energies of the maximally 
polarized state and the unpolarized state at / = 4/3, 
given by |P2/3-E'2 n |j wn l b e small compared to the energy 
difference due to the Zeeman field. Thus the total en- 
ergy difference between the fully polarized state and the 
pseudospin-balanced state is given by the Zeeman energy, 
which is cx B and favors the former, and the Coulomb 
energy, cx B 1 / 2 , which favors the latter. Consequently, 
there can be a transition between the two states, with 
the polarized state favored at larger magnetic fields. 

In neither of these states can there be skyrmions with 
negligble cost in anisotropy energy or Zeeman energy. 
Therefore, both states should have a reasonable energy 
gap, and should lead to well established quantized Hall 
states at low temperatures. 

C. / = 8/5 (i/ = -2/5) 

The situation at v = —2/5 should be very similar to 
that at v = —2/3. Again, we have two competing states: 
a balanced state where there are equal numbers of holes 
in K j" and K' \., and a maximally polarized state with all 
the holes in one of the two states. As there will be essen- 
tially no probability for two holes to be at the same point 
of space in either of these states, there will be no differ- 
ence in their valley-anisotropy energy. As at v = —2/3, 
the maximally polarized state will be favored by the Zee- 
man energy, while the balanced state should be favored 
by the Coulomb energy. Again, there can be a transition 



between the two states as one varies the magnetic field 
at fixed filling factor. 

D. Other fractions 

Quantized Hall states should be expected at many odd- 
denominator "Jain fractions" in the range —1/3 > v > 
—2/3 according to the composite fermion picture. In par- 
ticular, if one considers the integer QH state at / = 2 as 
the reference state, and attaches two Chern-Simons flux 
quanta to each hole inserted in that state, one expects 
to find quantized Hall states at fillings v = —p/ [2p ± 1), 
where p is an integer which describes the total number of 
holes per quantum of flux in the effective magnetic field 
B c ff = B/(l ± 2p). If p > 1, the holes can be distributed 
in various ways between the K f and K' 1 states. As we 
found for v = —2/5 and v = 2/3, which correspond to 
p = 2 with the two choices of ±1, we expect that states 
with different polarizations will differ very little in their 
anisotropy energies, but that states with greater imbal- 
ance will be favored by the Zeeman field whereas states 
with lesser imbalance will be favored by the Coulomb en- 
ergy. Thus we would expect a series of phase transition 
at all these filling fractions, as discussed in Ref. [2"2] . 

A value of v that is intermediate between two Jain 
fractions can be described by means of a fractional value 
of p. As discussed irP^l the ground state of such a sys- 
tem can be described in the composite fermion picture 
as state where an integer number of composite fermion 
LLs are filled for one of the two hyperspin species while 
there is a fractional occupation of the highest LL for the 
other species. When B is varied, phase transitions occur 
as the fractional occupation shifts from one species to 
the other. It is also worth noting that oscillations that 
may be interpreted as signatures of such phase transitions 
were observed^ even in region very close to v = —1/2 
where quantized Hall states were not apparent. 



E. Effects of spin canting 

In the previous discussion, we assumed that the FQH 
states with equal population of K,K' valleys are AF. 
Now we consider the effect of Zeeman interaction on such 
valley-balanced states, and show that it only leads to 
weak canting of spins in the direction parallel to magnetic 
field. This canting changes the energy of the state very 
slightly and does not influence the phase transitions be- 
tween valley-balanced and valley-imbalanced FQH states 
in the interval —2/3 < v < —1/3. In this Subsection, 
we also discuss possible candidate states at filling factors 
-Kv< -2/3. 

In order to consider the possibility of spin canting, in 
the regime 1 < / < 2, we must consider states in which 
the spins on the two sublattices can have orientations 
which are no longer opposite to each other. We shall 
assume the sublattice spins and the magnetic field direc- 
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tion to be coplanar, and we shall denote by 9x and 9x' 
the angles between the spins and the magnetic field, so 
that the angle between the spin orientations is 9x> — Ok- 
We define fx and fx 1 to be the filling factors on the two 
valleys, which, without loss of generality, we may restrict 
to satisfy < f K > < fx < 1, and f K + fx> = /• 

We wish to understand how various terms in the en- 
ergy depend on the the angles 9k and Ox 1 - For fixed 
values of Jk and fx 1 , the SU(4) Coulomb energy will 
be independent of these angles. The Zeeman energy, per 
flux quantum, will be given by 



valley-balanced state {Jk = fx 1 ) in the range 1 < / < 
4/3. 

If Jk 7^ fx' and if 7 is sufficiently small, one finds that 
the the lowest energy state will not have canted spins but 
rather will have collinear spins with 9a = and 9b = it 
(as we assumed in previous Subsection). In this case, 
spins of electrons in the K valley are aligned with the 
magnetic field, and spins in the K' valley point in the 
opposite direction, as we have assumed in our previous 
discussions. The sum of the Zeeman and anisotropy en- 
ergies for the collinear state will have the form 



£ z = -EzUk cos 9 k + fx 1 cos 9 K ') 



(18) 



To evaluate the energy due to the anisotropic valley in- 
teraction, we may make use of Eq. ( 16 ). The three terms 



in this equation should all be evaluated for the same an- 
gles 9k and 9k 1 ■ The term with total filling 2 — / should 
have sublattice fillings (1 — Ja) and (1 — Jb) while the 
term at / = 2 must have both hyperspin states filled. 
We have previously observed that the anisotropy energy 
for FQH states in the range < / < 2/3 should be ex- 
tremely small. Therefore, we may neglect the term with 
total filling 2 - / in Eq.([l6) and write, for 4/3 </ < 2, 



Ef{0 K ,0 K ,) = {f-l)E?(fi K ,0 K >) ■ (19) 

Since valley-anisotropic interaction is invariant under a 
uniform rotation of all electron spins, the anisotropy en- 
ergies can only depend on 9k> — 9k, the angle between 
the spins. Moreover, we may write"*' 



E: 



2 \?K' 



k) = -{^bV (9z + + coa{0 K > - K )]) 

(20) 

As we have assumed g± < 0, this energy is clearly mini- 
mized when 9k' — Ok — 7r, in which case we have 



Ef 



-(2^|)- 1 (/-l). 9z 



Ef 



(21) 



If fx = /if, and Ez is small but nonzero, the to- 
tal energy will be minimized by a choice Bk> = —Ok = 
(tt/2~ a) where the canting angle a is proportional to the 
quantity 7 defined in Eq.Q. More precisely, one finds 



cos 



X 



if 



4(/-l)' 



(22) 



provided that the right hand side is < 1. If the right hand 
side of (22) is greater than 1, one finds that Ok — 0k> — 



0; both spins are aligned with the magnetic field and the 
state is no longer canted. The sum of the Zeeman and 
anisotropy energies for the canted state may be written 

as 



£ an + E z 



Ef 



(tt) 



Ezf 



sin a . 



(23) 



so the decrease in energy due to the canting will be pro- 
portional to jEz, when a is small. Since 7 ~ 0.1, the 
energy gain from canting is indeed quite small for any 



E a 



£z = Ef(n)-E z (f K -f K ,). 



(24) 



Since the energy gain due to Zeeman term is linear in 
Ez, this state will have a lower energy than the canted 
state, for sufficiently small values of 7. We find that the 
collinear state, with 9k — and 9k> = 7r is stable if 



fx - fx > > 7 



fxf. 



if 



/-I 



(25) 



For a fixed value of / and 7, this inequality will be satis- 
fied provided that fx — fx 1 is larger than a critical value 
Sf c , which, for 7 = 0.1 and 4/3 < / < 2 will be well 
approximated by 5f c ~ 7/ 2 /(4/ — 4). For such states, 
the phase transitions between FQH states with fx = fie 
and fx 7^ /if will occur as a function of magnetic field, 
very similar to our discussion above. 

For fx - /if < Sf c , the ground state for fixed fx, fx> 
will have spins that are canted from each other, and the 
ground state energy will be slightly lower than it would be 
for the uncanted state. However the qualitative behavior 
of the energy as a function of fx — fx 1 will not be altered 
from that of ( 24 1 : the sum of the Zeeman and anisotropy 
energies will still decrease monotonically as fj± — fs is 
increased. Thus, in this case we would also expect to see 
the same series of transitions as a function of magnetic 
field, at fixed filling factor, with at most small changes 
in the locations of the transitions when one of the states 
involved has — fs < Sfc- 

We note that for 7 = 0.1, we find Sf c < 0.3 for all / 
in the range 4/3 < / < 2, and 5f c -4- 1 for / -> 2. Thus 
at / = 9/5 and / = 5/3, where one expects the ground 
state to have all holes in one valley, so that fx = 1 and 
/if = / — 1, we have Sf > Sf c , and the majority spins 
will be fully aligned with the magnetic field, as we have 
assumed in the previous sections. This will also be true 
for the unbalanced states at / = 12/7 and / = 8/5, where 
all holes sit in one valley. The balanced states at these 
filling fractions will be antiferromagnets with spins that 
are only slightly canted, which is again consistent with 
our previous assumptions. 

Now we turn to discussion of filling factor interval 1 < 
/ < 4/3. First, we note that the above picture might 
change completely if one considers filling factors where / 
is very close to 1. In this case the right-hand side of (22) 
may exceed unity, so for fx — fx 1 , the canted state will 
be replaced by a ferromagnetic state, where the spins are 
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fully aligned with the magnetic field, and the valley state 
is a singlet. For the spin-aligned state, one finds 

£ an + E z = Ef{ix) - 2(f - l)g^l 2 B )- 1 - fE z . (26) 

If 7 = 0.1, however, the state with fx = fx 1 will be fully 
polarized only for / — 1 < 1/39, and the canting angle 
will be relatively small {cos 9k < 1/2) for filling factors 
as close to 1 as / — 1 > 1/19. A state with maximum val- 
ley polarization, i.e., fx = 1 and fx< = f — 1 should have 
a fully aligned spins, for any value of / — 1. We have pre- 
viously noted that for typical experimental parameters, 
the lowest energy electron-like excitations should be skyr- 
mons with a small number of electrons in the minority 
valley and a small number of overturned spins. There- 
fore, it seems likely that the ground state at a very small 
positive value of /— 1 will be a state which is partially val- 
ley polarized and partially spin polarized. However, the 
ground state under these conditions would be a Wigner 
crystal rather than an FQH state. 

At filling factor / = 6/5, we expect that there should 
be a valley-balanced ground state with fx — fx' , whose 
spatial wave function would be the particle-hole conju- 
gate of the valley-singlet state at / = 4/5. However, we 
expect that the occupied hyperspin states at / = 6/5 
would be different, leading to a weakly canted AF state 
rather than a spin aligned state (as for f — 4/5). In 
principle, for large values of B, there could be a series of 
transitions to states with fx — fx' 0. However, it is 
likely that the Coulomb cost for such states will be high, 
and we would not expect to see such states for reasonable 
values of the magnetic field. At present no experimen- 
tal observation of an FQH state of any type has been 
reported at / = 6/5. 

We note that the parameter 7 will be increased if one 
applies a magnetic field parallel to the graphene sheet 
that is large compared to the perpendicular field. For 
sufficiently large values of 7 we expect that in all states 
with v < 0, the electrons in the zeroth Landau level 
will all have spins completely aligned with the magnetic 
field. In this case, one would expect that particle hole 
symmetry is restored between states / and 2 — /. 



VI. DISORDER EFFECTS 

The energy gap E g of a quantized Hall state may be 
defined as the energy cost to add a quasiparticle and a 
quasihole very far apart. In the limit of vanishing disor- 
der, the energy gap may be obtained from incompress- 
ibility or transport measurements. In transport measure- 
ments, the low temperature behavior of the longitudi- 
nal conductivity a xx is thermally activated, varying as 
e -E g /2 k B T _ Measurements of the compressibility using 
an SET- 20 * 22 -!, should show a discontinuity A/z in the elec- 
tronic chemical potential equal to E g e/e* at the quan- 
tized Hall density, where e* /e is the charge of the quasi- 
particle in units of the electron charge. (For all cases 



considered here, e/e* is equal to the denominator of the 
filling fraction /.) 

The effects of disorder on various types of measure- 
ments are only partially understood, even in the much 
studied case of GaAs systems. In principle, any amount 
of disorder should lead to a finite density of localized 
states in the energy gap, so that dfi/dn should be fi- 
nite even when n coincides with an ideal quantized Hall 
value, and the associated jump in /i should occur over a 
finite range of densities, inversely proportional to dfi/dn. 
Moreover, the peaks and valleys in /i on the two sides of 
the quantized Hall density, which would be sharp cusps 
in the ideal system, are typically rounded over a similar 
finite density range. In an SET measurement, the width 
of this range can be minimized, and the slope dfj,/dn, can 
be maximized, by placing the SET at a location that is 
relatively free of disorder. 

The behavior of dfi/dn described above can be seen 
in panel (a) of Fig.l below, which shows data for a sus- 
pended graphene sample near v = 0, at three values of 
B. Qualitatively similar behavior is seen in SET mea- 
surements near other quantized Hall states in graphene, 
as well as in GaAs. However, we do not have a good 
quantitive understanding of the density of states in the 
gap or of the extent of rounding in any of these cases. 

The effects of disorder on transport measurements are 
even more difficult to interpret. For an integer quantized 
Hall effect, in a system of non-interacting electrons, one 
might expect to see activated behavior at low tempera- 
tures, determined by a mobility gap, i.e., the energy dif- 
ference between extended states in the filled and empty 
LLs. However, in an interacting electron system, or in the 
presence of phonons, one would expect that the conduc- 
tivity will be described by variable range hopping at the 
lowest temperatures. Activated behavior is typically seen 
in experiments over an intermediate temperature range, 
typically about a factor of two in temperature, but the 
interpretation of the activation energy, and its relation to 
the energy gap of an ideal system are not at all clear 

Disorder effects in graphene may be particularly im- 
portant, and may present special features, in the cases 
where the low energy charged excitations are large-area 
valley skyrmions. Below, we shall explore the effects of 
disorder on quantum Hall states in graphene at different 
filling factors (both integer and fractional). 

The spin- and valley-polarized states that are predicted 
for various filling fractions in the range < / < 1 im- 
ply spontaneous breaking of a continuous SU(2) valley 
symmetry. As a result, these states will have low en- 
ergy Goldstone modes, and they will be very sensitive 
to valley-selective disorder potential. By the Imry-Ma 
argument!^, arbitrarily weak disorder will destroy the 
long-range valley order, and a ground state with spa- 
tially varying valley polarization should form. The pre- 
cise nature of the textured ground state, and the struc- 
ture of low-energy excitations depend on the dominant 
type of the valley-selective disorder. However, generally 
such textured ground states will have low-energy charged 
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excitations. This can significantly reduce the robustness 
of a quantum Hall state, or even completely destroy it. 

In contrast, valley-singlet states at < / < 1 (/ = 
2/5,2/3,4/9) do not have a spontaneously broken sym- 
metry. Therefore they should be fairly robust: weak 
valley-selective disorder might slightly reduce the value 
of the energy gap, and might introduce isolated localized 
states in the gap, but should not change the nature of 
the ground state. 

In the range 1 < / < 2, we find two types of quan- 
tized Hall states; both types turn out to be insensitive 
to weak disorder. First, states with unequal populations 
of the two valleys have a spontaneously broken valley 
symmetry, but the anisotropy energy reduces this from 
a continuous SU(2) symmetry to a discrete Z(2) sym- 
metry. The electron spins in the majority and minority 
valley states are aligned, respectively, parallel and an- 
tiparallel to the external magnetic field, so there will be 
no Goldstone mode associated with either the spin or val- 
ley degree of freedom. For states with a spontaneously 
broken Z(2) symmetry, the Imry-Ma argument suggests 
that disorder is marginal in two dimensions. Although 
domains should be induced, in principle, by any amount 
of disorder that couples to the difference in populations 
of valleys K and K' , the size of such domains should be 
exponentially large when the disorder is weak, so domain 
walls should be few and far aparPSI. 

Second, for states with equal population of valleys K 
and K' , we expect that the ground state will be a canted 
antiferromagnetic state, where the spins have a sponta- 
neously broken U(l) symmetry in the plane perpendic- 
ular to the magnetic field. This broken symmetry gives 
rise to a gapless mode in the spin sector, but the per- 
pendicular spin component would not couple to disorder 
produced by adatoms, impurities or strains. (In addition, 
even if one produces a skyrmion in antiferromagnetic spin 
order by forcing the direction of antiferromagnetic or- 
der to wrap around the unit sphere, such a skyrmion 
would not carry an electric charge.) The broken valley 



symmetry in an AF state is again pinned by the valley 
anisotropy terms, so that only a Z(2) symmetry is spon- 
taneously broken. Thus there are no Goldstone modes in 
the valley sector, so there would be no charged skyrmions 
with large area and low energy. Moreover, since the 
K and K' electron populations are equal, the broken 
Z(2) symmetry should not couple to potential fluctua- 
tions that break the valley symmetry. 

Because of this, all fractional states in the interval 
1 < / < 2, both valley-unpolarized and partially valley- 
polarized states, should be robust in the presence of weak 
disorder. For the case of valley-polarized states, spatially 
varying configurations of the order parameter can be cre- 
ated only if the disorder becomes stronger than valley 
anisotropy, and the valley-unpolarized states should be 
even more robust. 

We focus now on the valley polarized states with fill- 
ing fraction < / < 1, which are most susceptible to 
valley-selective disorder. It is convenient to specify the 
orientation of the valley polarization by a unit vector n(r) 
in pseudospin space, whose components describe the nor- 
malized expectation value of the operators ^ (r)fj^(r), 
for i = 1, 2, 3. We distinguish two types of valley-selective 
disorder: (i) an in-plane random field that couples to n x 
and/or n y , and therefore mixes valleys K,K', and (ii) a 
random field that distinguishes the two valleys and there- 
fore couples to n z . The first type of disorder could arise 
from adatoms situated at the center of hexagons or on 
hexagon sides in graphene lattice (e.g., various metallic 
adatoms). The second type of disorder may stem from 
adatoms that are situated on the sites of graphene lattice 
(e.g., H-adatoms), or from vacancies. Alternatively, such 
disorder can arise due to random short-range strain^} 
such strains are present in graphene on Si02 substrates, 
but probably are less important in suspended graphene. 
Due to the different origin of two disorder types, one or 
the other will likely prevail in any experimental situation. 



Assuming that the spin degree of freedom is frozen by the Zeeman interaction, we can write an effective energy 
functional for the SU(2) valley order parameter n(r): 

H =yJ dr I Vn ( r )] 2 + // drh ± (r)n ± (r) + fj drh z (r)n z (r) , (27) 

where nj^ = (n x ,n y ). We have introduced an explicit factor / in front of the random potential terms to emphasize 
that the coupling to the order parameter for a state that is fully valley polarized should be proportional to /, for a 
fixed value of B. , 



We shall assume that the random fields hi, h z , have 
short-range correlations on the scale of magnetic length. 
This is obvious for the case when disorder is due to 
adatoms, and it may be correct in many cases for lo- 
cal strains produced by inhomogeneities in the coupling 



to a substrate. Thus we may write 

(h a (r)) = 0, (h a (r)h a (r')) = V^n iniPi± S(r-r'), a = x,y. 

(28) 

(MO) = 0, (h z (r)h z {v')) = V*n imp , z S(r - r'), (29) 



where ni mPi2 and ni mPi j_ may be interpreted as densi- 
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ties of the two types of impurities, and Vq represents the 
root-mean-square value of the coupling to a single impu- 
rity. The quantity Vq has dimensions of an energy, and 
it will be proportional to the density of magnetic flux, 
£?/$()• Strictly speaking, Vq should be different for dif- 
ferent types of impurities, but we may ignore this distinc- 
tion without consequence for the following discussion. 

Following the Imry-Ma argument, the disorder poten- 
tial will destroy long-range valley order by creating do- 
mains. The domain configuration is determined by the 
disorder realization, and the structure of domains de- 
pends on the dominant type of disorder. The typical 
domain size £ can be estimated by noting that the en- 
ergy gained due to the adjustment of the order parameter 
to the local disorder configuration, E^ is ~ /Vo£^/Jii mp , 
should be comparable to the stiffness energy ~ p s . (No- 
tice that because of 577(2) symmetry, the order parame- 
ter varies on the scale £.) Thus, the typical domain size, 
given by 



i 



fV 0y /f, 



(30) 



is proportional to the valley-stiffness p s , and inversely 
proportional to the filling fraction. Here 7ii mp = 

max(jli m p jZ , TT-imp,_L ) ■ 

Let us consider two limiting cases: (1) n. imp _L 3> 
"■imp. 2, that is, in-plane disorder dominates, and (2) 
"imp._L "C Himp,2, the random h z field dominates. The 
two cases are characterized by different symmetry and 
different structure of the textured ground state, as dis- 
cussed below. However, in both cases there are low- 
energy excitations. 

To gain some intuition about case (1) (dominant in- 
plane disorder), let us first completely neglect h z . In this 
case, the order parameter in different domains will be 
pointing predominantly in-plane. This will lead to the 
formation of vortices, where the order parameter rotates 
by 2n. At the vortex cores, the order parameter must 
tilt out of the xy plane. Given the n z — > —n z symmetry 
of the Hamiltonian, there are two degenerate configura- 
tions of the order parameter for a single vortex (n z at 
the vortex core can be either positive or negative). Such 
configurations correspond to oppositely charged merons 
with charge ±e/2. This has two implications. First, the 
textured ground state has some density of frozen charges 
of each sign. Second, for an isolated vortex, changing 
the charge from positive to negative sign costs no stiff- 
ness energy. It should be noted, however, that the precise 
configuration of the order parameter will be determined 
by the balance between stiffness energy, disorder, and 
long-range Coulomb interaction. 

A weak random h z field will pin the direction of the 
n z component in the mcron cores in accordance with 
the local disorder configuration. As a result, the meron 
and anti-meron state will no longer be degenerate, but 
rather, will be separated by a small energy difference. 
The precise value of this energy difference will depend 
on the local disorder and on the Coulomb interaction 



with nearby frozen charges. Due to randomness, the en- 
ergy difference will be distributed in some energy interval 
(— E sk ; -E s fc); almost certainly, E sk will be much smaller 
than the skyrmion energy in a clean system, E® k = 4irp s . 
As we discuss below, this may strongly enhance the com- 
pressibility of certain fractional states, such that they will 
appear to be fragile or absent in experiments. 

Let us define a net skyrmion density iYsfe(r) as the 
difference in the number of negative and positive merons, 
per unit area. If we assume that the energy differences 
for converting a positively charged skyrmion to a negative 
one are uniformly distributed in the interval (—E sk ] E s k), 
we find that the energy cost to produce a net imbalance 
in the density of positive and negative merons should be 
given by 



E[Nsk\ = jj N 2 k (r)dr, 



(31) 



where K can be estimated in terms of the typical energy 
of a single extra skyrmion and the typical size of a single 
domain, 

K~E sk e. 

The chemical potential p for electrons is defined as the 
derivative of the total energy with respect to electron 
number, after subtracting any contribution from the in- 
teraction of the electrons with a macroscopic electrostatic 
potential. If we assume that fluctuations in the electron 
density arise purely from fluctuations in the skyrmion 
density, and are given by Sn = (e* /e)N sk , then we ob- 
tain an inverse compressibility given by 



dp 
dn 



(32) 



Taking into account Eq. ( 30 1 , we rewrite the above equa- 
tion as 



dp e 2 p 2 s E sk 
dn ~ VP 



i, "in.,. / 2 e* 



(33) 



Noting that E s k — Cp s , where C is a small numerical 
constant, we obtain the final form of inverse compress- 
ibility 



dp Ce 2 pi 

X 



dn V^n imp p e *< 



(34) 



This simple result allows us to make predictions re- 
garding the sensitivity of different ferromagnetic states 
at < / < 1 to disorder. For example, let us compare 
the integer state / = 1 (which is observed experimentally 
in compressibility measurements) to the fractional state 
/ = 1/3 (which is absent). The / = 1 is characterized by 
valley stiffness p s — e 2 /16V^Ib, while / = 1/3 stat e has 
a stiffness which is approximately 12 times smalleiP^ZI 
(note however that this number should be viewed only 
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as an estimate, due to large finite-size effects in numeri- 
cal calculations^^). Therefore, the ratio of the inverse 
compressibility for these two states is given by 



(dn/dn)\ f= i _ 12 3 
(dfi/dn)\f = i/ 3 81 



21. 



(35) 



Thus, in the disordered case the skyrmion contribution 
to compressibility is much smaller for the integer state 
/ = 1; this is due to the large valley stiffness of this 
state, which gives rise to large domains and smaller den- 
sity of frozen charges and low-energy excitations. This 
explains why the disorder-induced low-energy skyrmions 
may suppress the / = 1/3 state in experiment, but not 
the / = 1 state. 

Turning, now, to the case (2) above, where h z is the 
dominant random anisotropy term, we first note that ne- 
glecting in-plane disorder completely, the Hamiltonian 
(271 is characterized by a U(l) symmetry with respect 
to rotations of the order parameter in the xy plane. 
This symmetry reflects the fact that the energy does not 
change if we globally change the relative phase of K and 
K' states; such that if a(r)\K) + (3(r)\K') is a ground- 
state configuration for a given order parameter, then 
a(r)\K) + e lv f3(r)\ K') is also a ground state. An impor- 
tant difference compared to case (1) is that in the ground 
state there are no textures and charge excitations: the 
order parameter rotates only along some meridian of the 
Bloch sphere, not enclosing any area (such that topolog- 
ical index of the order parameter configuration is zero). 
However, the non-uniform configuration of the order pa- 
rameter in the ground state will be favorable for low- 
energy meron-anti-meron pairs, where phase <p rotates in 
different directions in neighboring domains. Therefore, 
low-energy excitations are expected in this limit as well. 

It is also worth noting that if small n- lmp ^ is added, 
larger super-domains will be formed in which the order 
parameter has a preferred component in the x-y plane. 
Such super-domains will enclose many domains with pre- 
ferred n z orientation. As a result, non-zero density of 
charges will be introduced in the ground state. It is 
not obvious whether the lowest-energy excitations will 
be textures in " small" domains, the size of which is set 
by riimp,z, or the merons/anti-merons forming in "large" 
domains, the size of which is set by nimp,±- I n either 
case, we expect that these excitations will have much 
lower energy than E° sk . 

The discussion above suggests that in the limiting cases 
when one kind of disorder dominates, there are low- 
energy charged excitations. In the intermediate regime, 
the textured ground state will have a certain density of 
charges (on the order of an electron charge per domain). 
Although it is difficult to make any quantitative predic- 
tion, it is reasonable to expect that the energy cost for 
creating skyrmions will be greatly reduced as well in this 
case. 

At this point, we are unable to say whether the en- 
hanced compressibility caused by skyrmions in the pres- 
ence of disorder is strong enough to account for the strong 



suppression of the FQH states with odd-numerator fill- 
ing fractions in the range < / < 1 in the SET mea- 
surements reported in Refs. [20122) . We also have little 
information about the mobility of these skyrmions, so 
it is difficult to estimate their contribution in transport 
measurements. Nevertheless, the apparent absence of an 
FQH signature at / = 1/3 in transport measurements is 
striking, and it seems likely that this absence is also due 
to the presence of low energy skyrmions at this value of 
/• 

It is also worth noting that suppression by disorder 
of quantum Hall states with a spontaneously broken 
sym metry has been observed in other multi-valley sys- 
tem j 46 | 47 | 49 | 50 l For example, in AlAs-2D electron sys- 
tem, which is characterized by a double valley degener- 
acy, the measured activation energy gaps of certain frac- 
tional states in valley-symmetric case can be reduced by a 
factor of 15 compared to the case when valley symmetry 
is broken by strairP^. It should be noted, however, that 
the quantum Hall ferromagnet in AlAs has a lower, Z-i 
symmetry compared to the quantum Hall ferromagnet in 
graphene at < / < 1, with the easy- valley anisotropy 
being generated by the asymmetry of the band struc- 
ture's. This may change the structure of domains, do- 
main walls and their excitations, compared to the case of 
graphene. We also note that experiments on Si (lll)p^ 
and Si (lOO^P^ have observed FQH states at various frac- 
tions, almost entirely with even numerators. However, a 
strong FHE state has been reported at nu = 1/3 in at 
least one experiment on Si (100)PM 



VII. EXPERIMENTAL MEASUREMENT OF 
ENERGY GAPS AT v = AND v = -1. 

Below, we present experimental measurements of the 
energy gaps of the broken-symmetry states at filling fac- 
tors v = and —1 in suspended graphene. We extract 
these energy gaps using a scanning SET. Details of the 
device fabrication and the measurement technique are 
provided in Ref. [20] ■ Briefly, we use a back gate to 
modulate the carrier density in the graphene, and mon- 
itor the resulting change in current through the SET to 
measure the local chemical potential \x. This is accom- 
plished using a DC feedback technique to maintain con- 
stant current through the SET; the change in sample 
voltage required for this feedback provides a direct mea- 
sure of the chemical potential of the graphene. 

Figure la shows the chemical potential as a function of 
carrier density at several representative magnetic fields. 
As the chemical potential approaches v = 0, it first de- 
creases, then abruptly rises, and finally declines again 
before stabilizing. The regions of decreasing chemical 
potential (negative compressibility) are understood^ to 
arise from electron-electron interactions, which become 
especially pronounced at the low electron and hole den- 
sities surrounding v — 0. Similar behavior occurs for 
v = — 1 (data not shown). The energy gaps at v = 
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FIG. 1: (a) Chemical potential \x as a function of carrier density n at magnetic fields B = 2.5, 5 and 10 T. The energy gap is 

given by the difference between minimum and maximum chemical potential, (b) Step in chemical A^i at filling factor v — as 

a function of magnetic field, (c) Step in chemical at filling factor v = — 1 as a function of magnetic field. The red curves in 

panels (b) and (c) are fits with B 1/2 scaling: Afj, = (23.1B 1/2 - 11. l) meV, for v = 0, and A/Lt = ( 8.5B 1/2 - 1.8) meV, for 
v = -1. 



and —1 arc defined as the difference between the mini- 
mum and maximum chemical potential around each in- 
compressible state, and they are plotted as a function of 
magnetic field in Figs, lb and c, respectively. 

The magnitudes of the energy gaps are significantly 
larger than has been previously reported^, which reflect 
both the high sample quality and also the low dielectric 
environment of suspended graphene. At B = 12 T the 
gap at v = — 1 is comparable to the theoretical skyrmion 
gap A„ = _i — ^ ss 23meV, estimated for dielectric con- 
stant e — 5.24. Energy gaps at both filling factors are 
well-modeled by B 1 / 2 scaling (red fits), consistent with 
expectations for interaction-driven states. Both fits, how- 
ever, require subtraction of a constant quantity, giving 
negative intercepts of a few meV. A negative intercept in 
the measured energy gap may arise in part from disorder 
broadening of the cusps in chemical potential, although 
we note that the v = offset of -11.1 meV implies sig- 



VIII. SUMMARY AND DISCUSSION 

In summary, motivate d by the puzzling results of re- 
cent experiment ;] 19 * 20 ^ , we studied fractional and inte- 
ger quantum Hall states in the zeroth LL of monolayer 
graphene. We argued that the behavior at filling factors 
— 2 < v < 1 and — 1 < v < is markedly different, 
owing to the different importance of short-ranged valley- 
anisotropic interactions. Qualitatively, one should think 
about fractional states at— 2 < v < — las electrons 
added to vacuum (empty Landau level). Since corre- 
lations in fractional states are such that the probabil- 
ity to find two electrons at the same point is negligible, 
the short-range valley-anisotropic terms do not play any 
role, which leads to an SU(2)-symmetric effective Hamil- 
tonian. In contrast, fractional states at — 1 < v < 
should be viewed as resulting from adding holes to the 



nificantly larger disorder than was estimated for other 
samples with comparable quality^. 

We would like to attribute the difference in energy gaps 
between v = and v = — 1 to the effects of skyrmions, 
which occur in the latter case. According to mean field 
theory, neglecting the effects of valley anisotropy and Zee- 
man energy, the energy gaps for creating a single electron 
or hole should be equal at these filling factors. However, 
as discussed in the text, the lowest energy charged exci- 
tations at v = — 1 should be large area valley skyrmions, 
which lead to an energy gap that is one-half the single- 
electron gap. At v = 0, we expect that large-area 
skyrmions will be suppressed due to the effects of val- 
ley anisotropy. The fact that the observed energy gaps 
at v = are more that twice the energy gaps at v = — 1 
may well be due to effects of residual valley-selective dis- 
order, which would affect the skyrmions at v = 1 most 
strongly, as discussed above. 

I 

v = state. Even though holes stay away from each 
other, they still interact with the background imposed 
by the v — state, which breaks SU(2) symmetry. 

The SU(2) valley symmetry at —2 < v < — 1 has 
two main consequences. First, fractional states with 
odd numerators must occur via spontaneous breaking 
of SU(2) symmetry. As a result, there are gapless 
neutral Goldstone modes that correspond to fluctua- 
tions of the valley polarization. The charge excita- 
tions are valley skyrmions, which have lower energy than 
the quasi-electron/quasi-hole excitations. Furthermore, 
the SU(2) symmetry of the Hamiltonian makes odd- 
numerator states very susceptible to weak disorder, which 
destroys the long-range valley order and introduces a 
certain density of skyrmions and anti-skyrmions in the 
ground state. In this case, the energy of adding another 
skyrmion/anti-skrymion is expected to be strongly re- 
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duced compared to the clean system. The lower energy 
skyrmion excitations and strong sensitivity to disorder 
make the odd- numerator states at— 2 < v < — 1 very 
fragile. 

Second, it is likely that at filling factors with even nu- 
merators in the interval —2 < v < — 1 ground states are 
5J7(2)-symmetric valley-singlets. In these states the low- 
est energy charged excitations are quasi-electrons/quasi- 
holes which can have either valley index, but are spin- 
aligned, and skyrmions, if they exist, would be energet- 
ically costly due to Zeeman energy. The valley-singlet 
states are not sensitive to weak valley symmetry-breaking 
disorder (since they cannot take advantage of adjusting 
valley polarization according to the local disorder con- 
figuration, unlike odd- numerator states). Thus, even- 
numerator states are expected to be quite robust. 

In the interval of filling factors —1 < v < 0, the valley 
and spin symmetries are broken by valley anisotropics 
and Zeeman interaction, which makes both odd- and 
even- numerator states robust, and leads to multiple 
phase transitions between states with different degree of 
spin and valley polarization. The v = state is a canted 
AF state, in which states \Ks), \K's'} with nearly oppo- 
site spins s,s' in the xy-plane, slightly canted in the z- 
direction, are occupied. The v = — 1/3 state is a partially 
valley- and spin-polarized Laughlin state of holes \K \) 
(or \K' \)) formed on top of a state with s =f, s' =1, 
which optimizes the combination of valley anisotropy and 
Zeeman energy. 

At v = —2/5 and —2/3 there are two competing states: 
a canted AF state with equal valley occupation (which 
can be thought of as a singlet state of holes on top of 
v = state), and a partially polarized state with un- 
equal valley and spin occupations. While anisotropy en- 
ergy is very close in these competing states, the polarized 
state can take advantage of the spin Zeeman energy by 
aligning spin direction of the majority component with 
the magnetic field. At low magnetic fields, Zeeman en- 
ergy is negligible, and the states with equal occupation 
of two valleys, favored by the Coulomb interactions, are 
expected to be ground states. At higher fields, when 
Zeeman interaction becomes increasingly important com- 
pared to the difference of Coulomb energies of the com- 
peting states, a phase transition into a partially polarized 
ground state is expected. Similarly, at other filling fac- 
tors at — 1 < v < 0, we expect phase transitions between 
composite fermion states with different degree of spin and 
valley polarization. 

By contrast, no phase transitions are expected at filling 
factors v = —4/3 or 8/5, because these states are valley- 
singlets and are spin-polarized along magnetic field such 
that they take maximum advantage of the Zeeman en- 
ergy. At filling factors v = —14/9,-11/7, and —7/5, 
phase transitions are in principle possible: at low B, 
minimally spin and valley-polarized states are favored 
by the Coulomb interactions (these correspond to filling 
lowest possible CF levels), while at high B maximally 
spin-polarized states will be favored by the Zeeman in- 



teraction. 

The above picture is consistent with the recent exper- 
iments^ 1 20 1 22 1 . First, odd-numerator states at — 2 < v < 
— 1 were eithe r com pletely absent or very weak in both 
cornpressibility^nHU and transport measurements^. This 
is likely due to effects of disorder which greatly increases 
compressibility. The robust even- numerator singlets have 
been seerl 2 ^, and no phase transitions as a function of field 
have been observed at this point. Second, robust od d and 
even- numerator states have been seen at — 1 < v < d 2 ° l 22 l, 
and a series of phase transitions as a function of magnetic 
field have been observed at all filling factors except for 
v = -1/3. 

Turning to integer states, we argued that at v = 0, the 
elementary charged excitations could be a single electron 
or hole, or could be a small skyrmion, which would have 
a similar energy cost. At v = —1, we argued that the 
negative charged excitation could be a large area valley 
skyrmion or a small area spin skyrmion, while the pos- 
itive (hole-like) excitation should be a large area valley 
skyrmion, which would reduce the size of the gap relative 
to that at v = 0. Effects of disorder should further re- 
duce the measured energy gap at v = — 1, relative to that 
at v = 0. This might explain the observed differences in 
the chemical potential jumps at these two filling factors, 
which are shown in data reported in this paper. 

Although we have assumed above that the parameters 
g± and g z fall in the range where v = state is in the 
AF phase, our principal conclusions will be unchanged 
if the ground state were the Kekule phase or the CDW 
phase. The main difference is that there would be no 
canted states in these cases. The states for < \v\ < 1 
would always have a majority spin parallel to B and a 
minority spin in the opposite direction, if the populations 
are unequal. 

Throughout the paper we have assumed the valley 
sensitive interactions to be instantaneous in time. For 
phonon-mediated interactions, this is correct to the ex- 
tent that the phonon frequency is large compared to the 
Coulomb exchange energy within a Landau level. Retar- 
dation effects could lead to corrections to our result that 
the energy of states with \i/\ < 1 should be independent 
of the orientation of the valley polarization < r >. 



IX. FUTURE DIRECTIONS 

Further insights into the nature of fractional states in 
the zeroth LL of graphene can be obtained from exper- 
iments with parallel magnetic field. At a fixed perpen- 
dicular magnetic field, it should be possible to observe 
phase transitions at various filling factors in the interval 
— 1 < v < as a function of parallel field component. 
Moreover, the dependence of the energy gaps on the par- 
allel field should reveal the spin structure of excitations 
at different filling factors. 

In the future, it would be also interesting to study the 
fractional quantum Hall effect in higher LLs of monolayer 
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graphene. Owing to the different structure of LL wave 
functions, the effective interactions will be different com- 
pared to the zeroth LL in graphene and to GaAs. This 
will change the phase diagram of fractional states, as well 
as their energy gaps. Valley anisotropies in higher LLs 
in monolayer graphene will be different from the zeroth 
LL, which may completely change the valley and spin 
structure of the fractional states. 

Furthermore, new effective interaction regimes are re- 
alized in related materials, such as bilayer graphene, 
which will almost certainly lead to interesting phenomena 
in the fractional quantum Hall effect regime. Theoretical 
studies suggest that band structure tunability in bilayer 
graphene may allow one to tune effective electron inter- 
actions by applying perpendicular electric fields. This 
can be used to stabilize the desired fractional states and 
induce phase transitions between them. Additional tun- 
ability can be achieved by changing the dielectric envi- 
ronment of the sample^. 

Finally, another promising direction is to study the 
structure of the FQH edge states in graphene. Owing to 
the atomically sharp confinement!^, it is possible to avoid 
edge reconstruction^!]. A recent theoretical study-23 sug- 
gests that this might enable the observation of universal 



Luttinger liquid behavior at the edges of FQH states. 
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